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, Abstract 

We perform a numerical analysis of the gravitational field of a global monopole coupled nonminimally to gravity, and 
find that, for some given nonminimal couplings (in contrast with the minimal coupling case), there is an attractive 
region where bound orbits exist. We exhibit the behavior of the frequency shifts that would be associated with 
O , "rotation curves" of stars in circular orbits in the spacetimes of such global monopoles. 
<D ■ 

7— 1 ; 

i PACS number(s): 11.27.+d, 04.40.-b, 98.62. Gq 

>; 

An interesting property of global monopoles is that they are configurations with energy density decreasing with 
the distance as r~ 2 |J. This is very suggestive in view of the fact that, naively at least, this is precisely what seems 
to be required to provide a natural explanation for the flatness of the rotation curves in spiral galaxies. Moreover by 
assuming the existence of a global monopole in a typical galaxy the total Newtonian mass contribution of the portion 
of the global monopole contained within r ga i is found to be M ~ ar ga i/2, where a — ^irGrj 2 is the deficit angle of 
the global monopole spacetime. If we take the radial extent of a galaxy to be r ga i ~ 15Kpc and consider a typical 
(j [ grand unified theory with a symmetry breaking scale r\ w 10 16 GeV, this mass turns out to be « 10 69 GeV which 
O h> is ten times the total mass due to the 10 11 solar-mass stars in the typical galaxy: M sta rs ~ 10 68 GeV. This again 
^ ' seems to be precisely what is called for in the observations [||. Finally, and what seems to be an added bonus, if 
we assume that the field of the monopole extends on average a distance of 10 galactic radii from the galaxy, where 
the configuration presumably coincides with that of the monopole centered in the neighboring galaxy, the total mass 
associated with each monopole (and thus with each galaxy) turns out to be 100 times that of the galaxy. Thus the 
total contribution of the monopoles to the average density in the universe is of the right order of magnitude to account 
for the inflationary prediction of a universe with critical density. This is in fact the basic argument that allows one 
to place bounds on the number density of global monopoles present in the universe || , the exact nature of which 
depends of course on the value taken for 77. Indeed global monopoles have been shown to be able to trigger inflation 
on their own for suitably high r\ 0, [H, Q. 

Unfortunately this picture doesn't hold up to scrutiny at the next level of analysis because of two problems: 1) 
upon substitution of Newtonian gravity by general relativity, it turns out [Q , that the linearly divergent mass has, 
at large distances, an effect analogous to that of a deficit solid angle plus that of a tiny mass at the origin. The mass 
of this "monopole core" is about 0.8 a. In fact Harari and Lousto M showed that this small core mass is negative and 
produces a repulsive potential. They studied the motion of test particles, in the large distance region, in the spacetime 
of a global monopole, concluding that there are no bound orbits; 2) the fact that the monopole configuration is rather 
unique, in the sense that it is basically independent of the ordinary matter content in the corresponding galaxy, which 
conflicts with the fact that there is a rather large range of galactic masses for which the dark matter component is 
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about 10 times more massive than the ordinary matter component ||. 

In spite of these remarks the initially mentioned results are in our view still very suggestive, specially if we look 
at the remarkable degree of universality that seems to emerge from the systematic study of galactic dynamics || : 1 ) 
the ordinary matter content is such a good indicator of the mass of the dark matter component; 2) the universality 
of the form of the radial distribution of the dark matter content; 3) other seemingly coincidental facts as scaling 
properties between dark and luminous galactic structure parameters. We, therefore, believe that the general idea of 
monopoles as seeds of galaxy formation, and models for the galactic (and perhaps cosmological) dark matter deserve 
to be further explored. 

The purpose of this Letter is to illustrate the possibilities in these models by considering the simplest nontrivial 
modification of the global monopole picture: nonminimally coupled global monopoles. In this work we will show that 
the situation regarding the problems alluded to above, can change dramatically with the introduction of a nonminimal 
coupling. 

Let's consider the simplest such possibility, namely a theory of a triplet of scalar fields <fi a , a — 1,2,3, nonmini- 
mally coupled (NMC) to gravity with global 0(3) symmetry which is broken spontaneously to U(l). The Lagrangian 
for the model is 



C = 



16tt 



R + F(R, cj> a cj> a ) 



(V0 a ) 2 + V((/) a (j) a ) 



(1) 



where V(<fi a (j> a ) is a potential that results in the breaking of the symmetry and is usually taken to be the Mexican hat 
potential V(<j> a (j> a ) = |(0 a ^ a - rff , and F is an arbitrary function of two variables. Note that in these theories the 
resulting potential for the scalar fields (the terms in the Lagrangian that are independent of the scalar field gradients) 
is 



(2) 



Therefore, the ordinary matter would affect the location of the minimum of the effective potential through its effect on 
R, thus avoiding the scenario where the monopole configuration is universal. This opens the possibility to recover the 
correlation between the masses in the dark and ordinary matter components of the galaxy that were alluded before. 

Next we focus on the existence or lack of existence of bound orbits in the corresponding spacetimes. To do this 
we will look at a specific example corresponding to the choice F — (£(f) a (j) a )R (where £ is a constant) and proceed to 
construct the resulting spacetime and analyze the motion of test particles. We will specialize to the case of a spherical 
and static solution to the coupled field equations. The configuration describing a monopole is taken as usual 



>> a = Vf(r) — , 
r 
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with x a x a = r 2 , so that we will actually have a monopole solution if / — > 1 at spatial infinity. We adopt the following 
metric: 



ds 2 = -N 2 {r)dt 2 + A 2 {r)dr 2 + r 2 d9 2 
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r sin 
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and analyze solutions of the gravitational and scalar fields equations describing global monopole configurations and 
the resulting space-time. Owing to the complexity of the resulting equations, we will perform a numerical analysis. 
To do so we adopt the following variables: 
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v{f) = \n[N{f)] , 



A(r) = 1 - a - 



2m(f) 



-1/2 



(5) 
(6) 



where 



1 + 2£A' 



A = 8-7TT7 2 



(7) 



where we have introduced the dimensionless quantity f = r/r c = r/A 
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The solutions are obtained by the standard shooting method S with boundary conditions ensuring a regular origin, 
i.e., /(0) = 0, m(0) = 0, m(0).f = —a/2, and that the spacetime is asymptotically flat but for a deficit angle jl(|, i.e., 
/(oo) = 1, i/(oo) = | In [1 — a]. The "shooting parameter" is taken to be f(0),f- 

Then we calculate the ADM (gravitational) mass of these configurations, (see JhJ for a rigorous definition of the 
ADM mass for spacetimes that are asymptotically flat but for a deficit angle) which can be easily evaluated from the 
integral 
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MadMq = lirn^oo m(r) = J ^47rf 2 T t *(f) - —J df 
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For all values of £ that we analyzed the resulting ADM mass was negative as in the minimal coupling case (here we 
must remind the reader that the ADM mass for quasi-asymptotically flat spacetimes does not share the positiveness 
property of the ADM mass of asymptotically flat spacetimes). For £ = — 2, for example (corresponding to the config- 
uration we will analyze in more detail) the resulting ADM mass is « —2.5 a. 

We now study the motion of test particles around a global monopole. To do this we consider the "equatorial 
geodesies" for the metric (|]) (i.e., geodesies lying on the plane 9 — n/2). We denote the 4- velocity of a test particle 
as = (t,r,0,(p), where the dot stands for derivation with respect to an affine parameter which in the case of a 
massive particle can be taken as the proper time. The resulting equation is 



exp[2z/(f )] 



(r) 2 + exp[2i/(f)] 



1 — a — 



2m(f) 



L 2 

(f) 2 



E 1 



(9) 



where E and L are the energy and the angular momentum per unit of mass of the test particle respectively. This 
equation shows that the radial motion is the same as that of a particle with a "position dependent mass" with energy 
E 2 /2 in ordinary one-dimensional, nonrelativistic mechanics moving in the effective potential 



Veff = exp[2i/(f)] 



L 2 
(f)2 



(10) 



In Figure [I] we show the effective potential V e ff (see eq. (10)) for the deficit angles a = 0.43 and a = 0.125 respectively 
and for three values of the angular momentum; we note the presence of a potential well and therefore the existence 
of stable circular orbits. 



The rotation curves (RC) of spiral galaxies are deduced from the red and blue shifts of the emitted radiation by 
stars moving in "circular orbits" on both sides of the central region |llj . We therefore evaluate the frequency shift for 
a light signal emitted from a "test star" in circular orbit in the spacetime of a global monopole and received by a static 
observer at infinity. The two maximum and minimum frequency shifts (associated with a receding or an approaching 
star respectively) turn out to be: 
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FIG. 1. Functional dependence of the effective potential V e // vs f for the case £ = —2, a = 0.43 (left fig.) and a = 0.125 
(right fig.). For each configuration we show three values of the angular momentum: L — 4 (solid line), L = 1 (dashed line) and 
L = 0.3 (dash-dotted line). Here r c = (ryA 172 ) -1 . 



z± = 1 



(1-a) 
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N(r)(N(r)-rN(r) t ?) 
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fiV(f) 
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with z defined as usual, 



Aui ujr — uje 



(12) 



where uje [ur) is the emitted (detected) frequency. We can define Zd = h(z+ — zJ) which would be the quantity 
that is operationally identified with the velocity |i>(f)| of the RC of spiral galaxies 
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N(f)(N(r) - rN(f),f). 
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In Figure || we show z+ (dashed line), z_ (solid line) and zjj (dash-dotted line) as functions off for the cases a = 0.43 
(left fig.) and a — 0.125 (right fig.) respectively. We note that even for this very simple model the figures that would 
correspond to the rotation curves contain a relatively "flat region" within the values of r corresponding to the stable 
orbits (i.e. the behavior of zd near its maximum). 

We are of course not suggesting that this simple model is in any way a realistic candidate for a natural explana- 
tion of the phenomena associated with the RC of spiral galaxies. However, the fact that the two main objections 
against the study of global monopoles in this context can be removed by such a simple modification as the intro- 
duction of a nonminimal coupling and the surprising coincidences mentioned at the beginning of this letter, provide 
very suggestive motivations for pursuing in the analysis of such possibilities. We must, nevertheless, emphasize that 
monopole-antimonopole pairs annihilate rather efficiently jl2j and therefore the viability of this scenario, as applied 
to the galactic rotation curves, will also depend strongly on the degree of asymmetry in the number of monopoles and 
antimonopoles in the early universe. 



In conclusion we have shown that the qualitative difficulties of the original monopole scenario for explaining the 
rotation curves of galaxies can be solved by the introduction of nonminimal coupling. The challenge for the future is 
to look for a specific realization of this more complicated scenario that combines the qualitative success of the latter 
model with the quantitative coincidences of the former. 
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FIG. 2. Functional dependence of z+ (dashed line), z_ (solid line) and zd (dash-dotted line) vs f for circular orbits in the 
case £ = —2, a = 0.43 (left fig.) and a = 0.125 (right fig.). The asterisk depicts the location of the radius beyond which the 
stable circular orbits cease to exist (r ~ 9r c ). 



ACKNOWLEDGMENTS 



We wish to thank R. M. Wald for helpful discussions and also to J. Guven for his comments. M.S. and D.S. 
would like to acknowledge partial support from DGAPA-UNAM project IN121298 and to thank the supercomputing 
department of DGSCA-UNAM. 



[I] M. Barriola and A. Vilenkin, Phys. Rev. Lett. 63, 341 (1989). 

[2] R. M. Wald (private communication). 

[3] W. A. Hiscock, Phys. Rev. Lett. 64, 344 (1990). 

[4] A. Vilenkin, Phys. Rev. Lett. 72, 3137 (1994). 

[5] N. Sakai, H. Shinkai, T. Tachizawa, and K. Maeda, Phys. Rev. D 53, 655 (1996); N Sakai, ibid. 54, 1548 (1996). 

[6] I. Cho and J. Guven, Phys. Rev. D 58, 063502 (1998). 

[7] D. Harari and C. Lousto, Phys. Rev. D 42, 2626 (1990). 

[8] M. Persic, P. Salucci, F. Stel, Mon. Not. R. Astron. Soc. 281, 27 (1996). 

[9] W. H. Press et al, Numerical Recipes: The Art of Scientific Computing (Cambridge Univ. Press, Cambridge 1985). 
[10] U. Nucamendi and D. Sudarsky, Classical Quantum Gravity 14, 1309 (1997). 

[11] V. C. Rubin, W. K. Ford, and N. Thonnard, Astrophys. J. 238, 471 (1980); V. C. Rubin, W. K. Ford, N. Thonnard, and 

D. Burstein, ibid. 261, 439 (1982); 289, 81 (1985). 
[12] D. P. Bennett and S. H. Rhie, Phys. Rev. Lett. 65, 1709 (1990). 



5 



